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SCATTERING  OF  PLANE  WAVES  FROM  A COMPLIANT  TUBE  GRATING 
IN  A LINEAR  VISCOELASTIC  LAYER 

INTRODUCTION 


The  use  of  compliant  tubes  as  a reflector  in  a fluid  medium  has  been 
studied  both  experimentally  and  analytically  by  Toulis,1  and  more  recently 
by  investigators  at  the  Naval  Underwater  Systems  Center  (NUSC).2’3  The 
latter  work  has  shown  that  through  the  use  of  multiple-layered,  compliant 
tube  baffles  of  wide  and  tunable  bandwidths  of  reflectivit>  can  be  achieved. 
Brigham  and  Radlinski4'6  have  developed  an  analytical  model  of  scattering 
plane  waves  from  an  infinitely  wide  array  of  elliptical  shells  in  a fluid 
medium.  The  calculations  from  this  model  agree  with  the  experimental  trans- 
mission measurements  from  both  plastic  and  steel  compliant  tube  arrays. 
Brigham7  has  extended  the  above  model  to  include  scattering  from  planar 
compliant  tube  arrays  in  an  infinite  elastic  medium. 

In  this  report,  the  scattering  of  a plane  wave  from  a planar  grating  of 
compliant  tubes  embedded  in  a linear  viscoelastic  layer  is  analyzed  through 
the  use  of  an  approximate  solution.  In  the  formulation,  the  tubes  are 
treated  as  rectangular  elements.  Only  the  case  where  the  tubes  are  aligned 
perpendicular  to  the  direction  of  the  plane  wave  that  is  normally  incident 
to  the  grating  is  considered  here.  The  reason  for  these  limitations  is  the 
severe  mathematical  complexity  of  analyzing  inclusions  of  arbitrary  shape 
and  orientation  in  a plane  layer.  Since  the  range  of  frequency  has  been 
extended  from  the  previous  works,  higher  order  structural  modes  are  important 
to  the  analysis.  The  solution  in  the  viscoelastic  layer  is  given  in  terms  of 
inhomogeneous  plane  waves  (that  is,  a wave  that  is  attenuated  along  the  lines 
of  constant  phase).  The  periodic  nature  of  the  solution  makes  it  possible  to 
reduce  a continuous  angular  spectrum  of  plane  waves  to  a discrete  number  of 
propagating  and  exponentially  decaying  waves. 

To  obtain  the  necessary  material  parameters  as  input  for  the  analytical 
models,  the  elastic  Young's  modulus  and  mechanical  loss  factors  were  measured 
in  the  100  Hz  to  2 kHz  frequency  range  using  the  Briiel  and  Kjaer  complex 
Young's  modulus  apparatus.  The  temperature-frequency  scaling  relationships 
were  applied  to  the  measured  data  to  calculate  the  elastic  modulus  in  the 
frequency  range  of  interest.8  Since  longitudinal  velocities  of  these  materials 
have  a relatively  small  variation  with  frequency  and  temperature,  measurements 
at  a single  temperature  and  frequency  as  reported  by  other  investigators9  were 
used  for  the  study. 
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Of  specific  interest  is  the  effect  of  the  elastomeric  stiffness  on  the 
compliant  tube  array  resonance  frequency  and  the  amount  of  absorption  for  a 
given  mechanical  loss  factor  of  the  elastomer.  The  dynamic  compliance  of  an 
individual  tube  in  the  array  is  also  a parameter  that  strongly  affects  the 
transmission  and  reflection  from  the  grating.  Comparisons  of  the  experimental 
measurements  and  the  predictions  of  the  model  are  given  for  broad  frequency 
spectra  and  wide  ranges  of  material  constants.  However,  they  are  not  meant 
to  include  all  measurements  made  at  this  time,  but  rather  to  illustrate  the 
capabilities  of  the  analytical  model. 


MATHEMATICAL  FORMULATION 
DESCRIPTION  OF  THE  PROBLEM 

This  study  concerns  the  scattering  of  a plane  wave  incident  through  a 
fluid  medium  on  to  a planar  grating  of  compliant  tubes  embedded  in  a visco- 
elastic layer.  In  particular,  insertion  loss  is  calculated  for  a plane 
wave  normally  incident  on  the  layer  in  a fluid  as  shown  in  figure  1. 


GRATING  IN  A LINEAR  VISCOELASTIC  MEDIUM 


In  a linear  viscoelastic  medium,  the  three-dimensional  wave  equation 
can  be  written  as 


(A+2p)'vA  -pyx  (2w)  = P 


(1) 


where  A and  p are  the  frequency  dependent,  complex  valued,  material  properties 
£hat  reduce  to  the  Lame/  constants  in  the  limiting  elastic  ca^e.  ^The^symbol 
u is  the  displacement  vector,  A = V • u is  the  dilatation,  2u  = V x u is  the 
rotation,  and  p is  the  density  of  the  viscoelastic  medium.  This  formulation 
is  equivalent  to  the  Kelvin-Voigt  model  with  harmonic  time  dependence  assumed. 
The  solution  to  (1)  can  be  found  by  using  a theorem  given  by  Helmholtz  which 
states  that  there  exists  a scalar  potential  <p  and  a vector  potential  A such 
that  the  displacement  field  may  be  expressed  as 


\T  = Vf  + Vx  A, 


where 


V.  A = 0 . 


(2a) 


(2b) 


2 


Figure  1.  Plane  Wave  Normally  Incident  on  an  Encapsulated  Grating 
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In  terms  of  the  scalar  and  vector  potential,  the  dilatation  and  rotation  can  be 
written  as 


A=  V • (V  *) 


2u  = v x V x A . 


Employing  the  usual  steady  state  conditions,  and  substituting  (2)  and  (3)  into 
Cl)  shows  that  the  equation  of  motion  (1)  is  satisfied  if 


V2  0 + k2  0=0 


V2  A + k2  A = 0 


where  the  complex  wavenumbers  are  given  by 


kp  =w/cp 


k,,  = w/c„  . 


Here  the  complex  velocities  are  defined  by 


cp  = P-v-  [(  * + 2p)/p]l/2 


c = P-v.  [p/p  ]1/2  . 


where  p.v.  denotes  the  principal  value  of  the  square  root  of  a complex  number. 


Since  the  resonant  compliant  tubes  have  lengths  much  greater  than  their 
transverse  dimensions,  the  waves  can  be  considered  to  be  essentially  two 
dimensional;  that  is,  the  wavefronts  are  perpendicular  to  the  xy  plane  and 
the  amplitudes  are  independent  of  z.  In  this  case. 


A = kAz  , 


where  k is  a unit  vector  in  the  z direction.  Consider  the  solution  to  the 
general  scalar  wave  equation  in  two  dimensions: 


(*L  +aL\ 
W V/ 


| + k | = 0, 


•P  IIUIII  I 
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where  k is  the  complex  wavenumber.  Separation  of  variables  implies  that  the 
general  solution  of  (7)  for  an  assumed  steady  state  condition  is  a general 
plane  wave  of  the  form 


iP*r  -Q*r 

e e 


where  P is  the  propagation  vector  and  Q is  the  attenuation  vector  or 

t t i(P  + iQ  )x  i(P  + iQ  )y  . 

$ =£  e v x ^x'  e y y)J 


Substitution  of  (8b)  into  (7)  yields 

k2  = (Px  + iQx)2  + (Py  + iQy)2- 


Now,  let 


k = (P  + iQ  ) 

X X X 

ky  ' (Py  + ‘V  • 


kx  = p.v.  (k2  - k 2)  1/2  . f9d j 

In  general,  the  components  of  the  propagation  and  attenuation  vectors 
must  be  determined  by  the  boundary  conditions.  In  the  cases  under  consideration 
where  the  compliant  tubes  are  arranged  in  a planar  grating,  which  for  computa- 
tional purposes  is  assumed  infinite  in  extent  along  the  y-axis,  these  vector 
directions  are  readily  determined. 

If  the  solution  is  periodic  in  y,  then  (for  normally  incident  waves) 


£(x,y)  = I (x,  y+d), 

where  d is  the  periodic  grating  spacing.  Substitution  of 
to  the  condition 

iP  d -Q  d 

e ye  y = l. 


(8b)  into  (10)  leads 


(11a) 


From  (12),  one  sees  that  the  attenuation  is  always  in  a direction  normal 
to  the  plane  of  the  grating.  When 

Rek2>(^)2.  (Ua) 


the  waves  are  propagating  waves  whose  direction  is  determined  by  the  order  n. 
For  n equal  to  zero,  the  direction  of  propagation  and  attenuation  is  the  same. 
For  other  values  of  n,  the  amplitude  of  the  wave  varies  along  the  wavefront. 
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are  referred  to  in  the  elastic  or  fluid  medium  as  the  evanescent  modes  for 
which  the  waves  propagate  along  the  grating  and  decay  exponentially  normal 
to  the  grating.  As  seen  in  (12),  the  attenuation  of  evanescent  waves  in  a 
linear  viscoelastic  medium  is  normal  to  the  grating;  thus,  the  propagation 
necessarily  cannot  be  normal  to  the  attenuation  and  is  the  only  wave  that 
cannot  propagate.  For  more  complete  treatises  on  linear  viscoelastic  media, 
consult  the  works  of  Cooper,10  Cooper  and  Reiss,11  Buchen,12  and  Borcherdt . 1 3 ’ 14 

Assuming  plain  strain  conditions,  the  normal  and  tangential  deflections 
are  calculated  from  scalar  and  vector  potentials  as 


u = Vp  + V 

x A ; 

(14a) 

then 

u =££.  + 
x ax 

dA 

z 

9y 

(14b) 

u « !£.  - 
y ay 

dA 

' 

(14c) 

The  stress  is  calculated  from 

cr ..  = AA<5  ..  + H-  ( 
iJ  xj  \ 

' 8u. 

8u.  \ 

8xi  / 

(15) 

where  the  Kronecker  delta  6 is  defined  as 

6ij=1 

i = j 

6ij=° 

i ^ j, 

and  the  dilatation  A is  given  by 

dux 

A - A 

au 

a x 

+ ~dj 

ANALYTICAL  DESCRIPTION  OF  COMPLIANT  TUBES 

Usually  described  as  a long  oval  shell,  deflection  profiles  for  bending 
and  membrane  modes  of  a compliant  tube  are  exemplified  in  figure  2.  As  shown 
in  an  analysis  by  Toulis1  and  confirmed  at  NUSC,3  the  frequency  and  waveform 
of  the  fundamental  even  bending  mode  of  a compliant  tube  can  be  estimated  by 


JsMi. 
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BENDING  MODES 


FUNDAMENTAL 


SECOND 


FIRST  HALF-CYCLE 
SECOND  HALF-CYCLE 
UNSTRAINED  OVAL  CENTERLINE 


MEMBRANE  (HOOP)  MODE 


Figure  2.  Deflection  Shape  of  the  Lowest  Membrane  and  Two  Lowest 
Bending  Modes  of  a Compliant  Shell 


8 


TR  5433 


In  the  present  formulation,  the  compliant  tube  is  treated  acoustically 
as  long  thin  rectangular  plates  of  width  2b  with  rigid  connections  of  width 
2a,  but  is  modified  structurally  to  account  for  the  decreased  compliance  with 
decreased  eccentricity.  Since  only  plane  waves  normally  incident  to  the 
grating  are  considered  here,  only  modes  symmetric  about  the  x-axis  are  excited. 

At  the  surface  of  the  compliant  tube,  normal  stress,  tangential  shearing 
stress,  and  radial  and  tangential  displacements  are  continuous  for  bonded 
elastomers.  For  thin-walled  tubes  and  small  displacements,  the  stress  boundary 
conditions  are  equivalent  to  radial  and  tangential  force  equations  defined  in 
terms  of  the  shell  displacements.7  All  stresses  are  taken  at  the  mid-surface 
of  the  shell.  For  certain  conditions,  the  shell  equations  for  the  compliant 
tube  can  be  decoupled.  At  low  frequencies  where  only  bending  is  important, 
the  inextensional  condition  can  be  used  to  derive  a single  differential  force 
equation  relating  the  normal  deflections  to  the  external  stresses.  To  obtain 
the  membrane  inodes,  a single  differential  equation  results  by  assuming  that 
the  motion  is  primarily  radial  and  thus  setting  tangential  displacement  to  zero 
and  neglecting  bending  terms.  Furthermore,  for  thin  walled  cylinders,  the 
effect  of  the  tangential  shearing  stress  is  much  smaller  than  the  normal  stress. 
A simplified  equation  for  the  normal  bending  is  given  by  the  plate  equation 


ei>2  d“v 

12(1- dy4 


(16) 


where  a is  the  applied  normal  stress,  V is  the  normal  displacement,  E is  the 
Young' s modulus  of  the  shell  material,  h is  the  shell  thickness,,  ps  is  the 
shell  density,  and  v is  Poisson's  ratio.  The  rigid  body  motion  is  found  by 
simply  setting  E to  zero.  The  above  equation  for  bending  is  further  simplified 
by  use  of  the  normal  modes  defined  for  the  homogeneous  form  of  equation  (16) 
such  that 


where  ip  is  the  eigenfunction  of  the  normal  bending  mode  with  the  corresponding 
amplitude  V . u>t  is  the  resonance  frequency  of  the  t-th  normal  mode.  By  using 
(16),  a small  correction  to  the  effective  inertia  of  the  shell  wall  has  been 
neglected  for  eccentricities  not  equal  to  one.7 


The  fundamental  bending  mode  reflection  profile  for  a clamped-clamped 
beam  of  length  2b  is  taken  as 


= NeJl-^)  -b<y<b. 


(18) 


Here,  Ne2 


is  the  normalization  constant  calculated  from  the  condition 


,2 

4> 2 1 dy  = b, 


(19) 


which  results  in  Ne2  = 1.11.  For  compliant  tube  eccentricities  that  are  small 
compared  with  1,  the  deflection  profiles  of  an  elliptical  ring  as  given  by 
Brigham15  are  used  for  ^ . The  deflection  profiles  for  the  fundamental  even- 
even  bending  mode  for  various  values  of  ratio  a/b  are  shown  in  figure  3.  As 
a/b  approaches  one,  the  net  deflection  approaches  zero,  and  the  tubes  become 
less  compliant.  Either  a Fourier  series  fit  or  least  squares  fit  to  the 
deflection  profile  is  used.  Higher  order  bending  modes  that  are  resonant 
in  the  frequency  range  are  also  included  in  the  analysis.  The  rigid  body 
wave  function  to  a first  approximation  is  a constant  which  is  normalized  as 
above . 

For  plastic  compliant  tubes,  the  fundamental  membrane  mode  is  important 
to  the  analysis,  since  it  becomes  resonant  in  the  frequency  ra’ige  of  interest. 
The  eigenvalues  and  eigenfunctions  of  the  fundamental  membrane  mode  of 
uniformly  thick  elliptical  shells  have  also  been  determined  by  Brigham.16  The 
resonance  frequency  of  the  fundamental  membrane  mode  is  given  by 


7TCS 
wi~  2S 

o 


nfe/g  (1-v2)]  1/2 

2S 

o 


(20) 


Here,  c is  the  speed  of  sound  in  the  compliant  tube  material  and  Sq  is 
defined  by  tfie  elliptic  integral 


So'b 


/’ 

•'A 


[1-e 


2sin2  0 


.1/2 


d0 


(21) 


where  e is  the  eccentricity 
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If  the  tangential  displacement  is  assumed  zero  and  bending  terms  assumed 
negligible  near  the  membrane  resonances,  then  the  resulting  differential 
equation  describing  membrane  motion  is  given  by  the  reduced  radial  shell 
equation 

r Ehs  . i 

lr2,l-v2,  SSJ  (22) 

where  r is  the  local  radius  of  curvature  of  the  shell.  Again  using  the  normal 
modes  defined  for  the  homogeneous  form  of  (22)  such  that 


’ Eh 

s 

2 , 2 
Lr  (1-v  )J 


2 V i p = p h 2 V w & , 

m*m  Hs  s m w m 


we  then  have 


P h Y' ~ w2)  V i = ,t 
s s ^ m ’ m °xx  . 


DESCRIPTION  OF  COMPLIANT  TUBE 
GRATING  IN  A THIN  VISCOELASTIC 
LAYER  IMMERSED  IN  A FLUID 

The  various  regions  in  which  a solution  for  a normally  incident  plane 
wave  can  be  found  by  matching  the  appropriately  chosen  boundary  conditions  at 
the  interfaces  are  shown  in  figure  1.  The  solutions  are  chosen  to  satisfy 
the  condition  that  U is  zero  at  the  planes  of  symmetry  shown  in  figure  1. 

The  solution  for  theypressure  in  the  field  (region  1)  is  given  as 


n - jk  x 
Pj  - e o 


^ r -iki^x 


(24a) 


P = T P e n 
r *-•  n 


ffl- 


(24b) 


where 


(24c) 
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Here  is  the  incident  plane  wave  pressure  and  P is  the  reflected  wave. 
The  wavenumber  is  kQ  in  the  fluid  and  d is  the  periodicity  of  the  grating. 

The  viscoelastic  layer  is  divided  into  three  regions,  as  shown  in 
figure  1.  The  scalar  and  vector  potentials  for  the  dilatational  and  shear 
waves,  respectively,  are  given  as 


Region  2 


ik<2>* 
“P 


e 


(25a) 


where 


lk<2>x 


s 


(25b) 


• wT 

>?■[>.■  -W]“ 


Here  k is  the  complex  dilatational  wavenumber  and  k is  the  complex  shear  wave- 
number?  s 


Region  3 

The  potentials  in  this  region  are  chosen  to  satisfy  the  additional  con- 
straint that  U is  zero  at  y = b.  In  other  words,  for  a normally  incident 
dilatational  wave,  there  is  no  translational  motion  of  the  compliant  tubes 
in  the  y direction: 


0 = 
3 


E K 

n=0  l n 


(3) 


ik(3)x 
e n 


n p 


| cos  2 n 7T  [(y  - b)/|  j (26a) 


A3 


Y lC  (3)„ik!3)x 

n=0  n 


(3) 


e n 


+ Dn(3)eiknsX  | sin2n7r^(y  - b)/|  j, 


TR  5433 


where  t = d-2b  and,  to  satisfy  the  wave  equation. 


Region  4 


k <3>  = k 2 

np  lP 


k <3>=k2 
ns 


-(t): 

- w 


2l  1/2 


2 1 1/2 


(26b) 


(27a) 


where 


d'-’-wr 


(27b) 


Thus,  the  solutions  in  the  viscoelastic  layer  and  fluid  are  taken  as  the  sum 
of  propagating  and  evanescent  plane  waves  whose  direction  is  determined  by 
the  periodicity  of  the  layer.  Of  course,  the  higher  the  order  n of  an  evan- 
escent wave,  the  more  rapid  the  decay  away  from  the  grating. 

The  conditions  in  region  5 can  be  chosen  as  another  fluid  medium,  an 
infinite  impedance  medium,  or  a finite  elastic  plate  backed  by  a light  or 
heavy  fluid  medium.  The  chosen  backing  then  determines  the  boundary  condi- 
tions at  the  interface  between  regions  4 and  5. 
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In  the  fluid  region  behind  the  layer,  designated  as  region  5 in  figure  1, 
the  pressure  is  taken  as  a sum  of  outgoing  plane  waves  with  the  periodicity 
of  the  grating,  that  is. 


v £p 

• — A 


£5>  » 


t 

’ e p cos 


(28a) 


where,  again 


1/2  • 


(28b) 


BOUNDARY  CONDITIONS  AND  INFINITE 
SYSTEM  OF  EQUATIONS 

The  boundary  conditions  to  be  satisfied  at  the  interfaces  between  the 
designated  regions  are  as  follows: 

A.  At  x - -t. 


1.  The  tangential  shear  stress  in  the  elastomer  must  be  zero: 

rr(2)-n 

°xy  = °* 


which  in  terms  of  the  potentials  gives 

12) 


(29a) 


(29b) 


2.  The  normal  stress  a in  the  elastomer  must  be  equal  to  the  negative 


of  the  total  pressure  in  the  fluid 


£>  - -pt  - -(P( . pr, 


(30a) 
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or 


n=o  p J 


^(cfe-^i-Df  .“C4)}- (a^) 

-ik  t.  y*  ik^t,  /_  \ 

o 1 , nr  n 1 1 2nn  \\ 

+ n=0  Pn  6 cos  \d)  * 


4 ir  n 
d 


= e 


(30b) 


3.  The  normal  deflection  at  the  viscoelastic-fluid  interface  must  be 
continuous , 


.(2) 


dP_ 


o)2p  dx 


(31a) 


or 


.IK’ 

-£  ikm  Pj  e‘k™  ‘l 


(3lb) 


B.  At  x = x.  = +a 
J 


1.  The  normal  displacements  of  the  viscoelastic  layer  are  continuous 
across  the  elastomer  and  equal  to  the  compliant  tube  displacement  at  the 
tube  walls: 


■ 
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which  results  in 


t { (a«>  y *. . b«>  ey^y^>  y’ 

+ B«>  e-^^cos  (-27) 

- £ { <>  (a<?>  e<3pS  -B«  e-<3pS) 


/ 2nn\  /1(3)  ik^x.  , p.  (3)  -ik^x.\\  /y-b\ 

I — — I IC'  ' e n j + Dwe  n i|>  cos  2niT  I- — I 

\l  J\n  s 1 n s }Jf  J 


where  x„  = 


-a  and  = a,  and 


u(i,=Z.V®<;  y<|b| 

X t t t 


leads  to 


(32b) 


(32c) 


SK(\®  'lk"pX'-Bnfl)  xi) 

+ (t)(c»U)  * V1’ ^»!S )}«(*■)  - f v®*t 

y< |b| . 1 

The  amplitudes  of  the  tube  displacement  differ  at  most  by  a sign. 

2.  The  equation  of  motion  of  the  compliant  tube  is  satisfied  by  the 
approximation 


(32d) 


2 p h 2- 

s s 


W ~) 


v» 


, <2>. » <4> 

XX  - XX 


y<  lbl 


(33a) 
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2«A  Jj  (",2-“'2)vt 

=-|0{[X  ) 2](An’  Bf  •*;•) 

+ 1^1  *»  [cf  e-‘k!f’a  - Df  e*”*  ] }coS 

4^  -£>  ( <*>.*£*  - B<V<48’a  ) }c„8  (^f  ) 


*JsHxk» + 2(1 


where  t includes  all  modes  symmetric  to  the  x-axis. 


(33b) 


For  even-even  inodes  the  motion  of  the  tube  walls  is  symmetric  with 
respect  to  the  y-axis.  For  even-odd  modes,  the  motion  of  the  tube  wall  is 
antisymmetric  with  respect  to  the  y-axis. 


The  normal  stress  is  continuous  across  the  elastomer: 


0)  = (3) 

axx  ®xx 


b <y  j = 2,4 


(33c) 
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oOV^oosfe)  - £{[-2P 

J ' ' n=0  L ' p 


+ b'V^SI  -if  «<3’  [v‘k"?Xi-Vik"l>Xi  ]}oos2n,[(y-b,/l].  (33d) 
n P J * s L 

3 Tangential  displacements  are  continuous  across  the  elastomer  and  for 
bending  mode?  the  displacements  are  equal  to  the  inextensional  tangential  dis- 
placement across  the  tubes: 


yields 


u*  - U(3> 

y y 


b <y< 


d/2  j = 2,4 


i (\u'elkv  + 

n=0  ' ' 

t { 2ffn,3,-p’Xi  +B<3»a-<3pS) 


(34a) 


(34b) 


and 


(j)  = 0V^_  = r J_sv(j)i|/t  y<|b| 

y 8y  8y  t t 


(34c) 
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gives 


( 2ira  /.(])  ik^3  x.  -,0)  -ik^]  x.  \ 
< - —T  l A e n j + B e n j I 
\ d y n pJ  n P J / 


= rZV, 


(j)  d^t 


t 3y 


(34d) 


For  membrane  modes,  the  tangential  displacement  of  the  compliant  tube  is 
assumed  to  be  zero. 


4.  The  tangential  shear  stress  is  continuous  across  the  elastomer: 


axy(j)  " °xy^  |bky<|d/2|  3 = 2,4 


(35a) 
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and 


(i)  _ 


xy 


= o 


y<lb 


(35c) 


\ 


gives 


«|0K;  (1a£)[vlk"pXj  - 

i(k% ) 2 - W 1 [c<» + ^ ] h ( t)  - «• 

For  thin  walled  compliant  tubes,  it  is  assumed  that  the  tangential  shearing 
stress  does  not  produce  a comparable  stress  in  the  tube  walls. 

C.  At  the  elastic-fluid  interface  (x  = t2) 

1.  The  tangential  shear  stress  is  zero: 


«r(4)  = 0 * 

xy 


(36a) 


or  in  terms  of  the  potentials 


(<4>  - <4)  e-lki» 

+ k:  )2  - W ] ^ + D^ik"4*>t3) } - m - • 

2.  The  normal  stress  in  the  plate  is  equated  to  the  negative  of  the 
transmitted  pressure  in  the  fluid: 


(36b) 


„ <4) 

Pt  = axx 


(37a) 
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gives 


- t ik  <5)t„ 
y P e n 3 cos 

n=0  n 
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3.  The  normal  deflection  at  the  elastic  interface  is  continuous 


(37b) 


1 9Pt 
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(38b) 


The  above  equations  are  multiplied  by,  say,  the  orthogonal  functions 
sin  (2imr/d)y  or  cos  (2miT/d)y,  and  integrated  over  the  typical  grating  spacing 
of  0 to  d/2.  The  equations  defined  by  (33a)  are  multiplied  through  by  the 
appropriate  deflection  profile  <Pn  and  integrated  about  the  tube  in  the  interval 
y £ b . For  the  comparisons  of  experimental  data  with  calculations  as  shown 
in  the  next  section,  the  analysis  has  been  simplified.  First,  the  tubes  have 
been  considered  as  infinitesimally  thin  (that  is,  a-K)) , and  thus  the  width 
and  contribution  of  region  3 shrinks  to  zero.  Second,  since  only  tightly 
packed  compliant  tube  arrays  are  considered,  a good  approximate  solution  of 
closed  form  is  obtained  by  satisfying  the  displacement  equations  over  the 
entire  grating  spacing  d at  x = 0 and  the  stress  equations  only  over  the  tube 
walls.17  As  the  grating  spacing  increases,  full  consideration  must  be  given 
to  the  interstices  between  the  tubes. 
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EXPERIMENTAL  PROCEDURE 

Comparisons  of  transmission  through  the  layer  of  compliant  tubes  in  a 
viscoelastic  material  were  made  with  both  0.375  in.  (9.525  mm)  major  axis 
Lexan  plastic  tubes  and  0.925  in.  (23.5  mm)  major  axis  steel  tubes.  The 
plastic  tubes  had  an  a/b  ratio  of  about  0.5  and  the  in-air  fundamental  flexural 
resonance  was  measured  to  be  approximately  22  kHz.  The  static  compressibility 
of  the  plastic  tubes  was  20  times  greater  than  thai  of  water.  A cutaway 
plastic  tube  panel  is  shown  in  figure  4.  The  grating  spacing  in  this  example 
is  about  0.45  in.  (11  mm)  and  the  panel  is  of  dimensions  27  in.  x 27  in.  x 
0.375  in.  (0.69  m x 0.69  m x 9.525  mm).  The  14  in.  (0.35  m)  long  steel  tubes 
had  an  a/b  ratio  of  about  0.1,  an  in-air  bending  resonance  of  7.0  kHz,  and  a 
static  compressibility  of  about  one  hundred  times  of  water.  The  adhesive  bond 
between  the  steel  tubes  and  the  encapsulants  was  found  to  be  stronger  than  that 
with  the  plastic  tubes. 

Measurements  were  initially  made  of  the  transmission  through  the  grating. 

A standard  receiving  hydrophone  was  placed  in  close  proximity  to  the  center  of 
the  array  on  the  panel's  transmission  side  So  measure  the  insertion  loss. 
Insertion  loss  as  used  is  defined  as  20  log^ j Transmitted  Pressure/Incident 
Pressure | . Insertion  loss  measurements  were  made  to  determine  the  effect  of 
the  elastomer  on  shifting  the  array  resonance  frequencies  and  decreasing  the 
Q value  (quality  factor)  of  the  resonances.  The  specific  impedance  pc  of  the 
encapsulants  in  all  cases  was  close  to  that  of  water.  Here  c is  the  longitud- 
inal velocity  of  the  material . 

INSERTION  LOSS 

The  effects  of  the  viscoelastic  encapsulant  on  the  array  resonances  are 
considered  in  this  section.  In  figure  5,  a comparison  of  the  measured  and 
calculated  insertion  loss  is  made  for  the  plastic  tubes  embedded  in  a low 
stiffness,  low  damping  castable  rubber,  manufactured  by  B.  F.  Goodrich,  that 
has  both  a density  and  longitudinal  velocity  close  to  that  of  water.  The  0.5  in. 
(12.70  mm)  thick  panel  had  a 0.45  in.  (11.50  mm)  grating  spacing.  The  funda- 
mental bending  mode  and  the  first  membrane  mode  are  the  most  significant  modes 
in  the  measured  frequency  range. 

The  shear  modulus  of  the  Goodrich  material  was  determined  to  be  200  psi 
(1,38  x 106  n/m2)  at  10  kHz  with  a corresponding  loss  factor  of  0.28.  The.' 
loss  factor  is  the  ratio  of  the  imaginary  part  of  the  complex  modulus  to  the 
real  part.  These  values  were  essentially  constant  in  the  temperature  arid 
frequency  ranges  of  interest . Close  agreement  is  found  between  the  measured 
and  calculated  insertion  loss  in  figure  5.  The  measured  fundamental  array 
bending  mode  resonance  at  about  18  kHz  is  similar  to  that  of  a single  plastic 
tube  in  water.  The  measured  insertion  loss  does  not  vary  significantly  from 
that  of  an  unencapsulated  array  of  similarly  spaced  compliant  tubes  found  in 
reference  6,  where  at  the  array  resonance,  the  unencapsulated  panel  acts  as  a 
dynamically  soft  reflector. 


23 


• -v 


TR  5433 


The  plastic  tubes  were  then  encapsulated  in  a 0.375  in.  (9.5  nun)  thick 
polyurethane  panel.  Measurements  of  the  moduli  of  the  polyurethane  PRC-1534 
could  not  be  made  because  it  contained  a carcinogen  and  is  no  longer  com- 
mercially available.  A typical  value  of  the  shear  modulus  of  the  polyurethane 
PRC-1564  was  determined  to  be  about  20,000  psi  (1.38  x 10s  n/m2)  with  a cor- 
responding loss  factor  of  0.48.  The  shear  modulus  of  polyurethane  is,  however, 
strongly  temperature  dependent.  The  calculated  and  measured  insertion  loss 
curves  for  the  polyurethane  panel  at  22°C  are  shown  in  figure  6.  The  closely 
spaced  tube  array  encapsulated  in  this  higher  shear  modulus  material  provides 
a better  impedance  match  with  water  and  thus  more  energy  is  transmitted  through 
the  panel.  Also,  the  higher  shear  modulus  contributes  to  stiffness  loading 
which  increases  the  frequency  of  the  array  resonance.  Because  of  its  low  shear 
modulus,  the  attenuation  per  unit  length  of  the  Goodrich  material  is  greater 
than  that  of  the  polyurethanes. 

Comparison  of  the  measured  and  calculated  insertion  loss  for  the  steel 
compliant  tubes  with  a 1 in.  (25.4  mm)  grating  spacing  and  encapsulated  in 
approximately  a 1 in.  (25.4  mm)  thickness  of  the  low  shear  rubber  is  shown 
in  figure  7.  In  this  example,  the  array  radiation  loading  of  the  closely 
spaced  tubes  alone,  results  in  the  array  resonance  being  at  a higher  frequency 
than  that  of  a single  tube  in  water.  The  resonance  of  the  first  harmonic  even 
bending  mode  is  seen  at  about  40  kHz  in  both  the  measurements  and  calculations. 
These  curves  are  to  be  compared  with  the  insertion  loss  for  the  same  array 
encapsulated  in  NUSC-A  polyurethane  as  shown  in  figure  8.  The  polyurethane  is 
a castable  mixture  of  100  parts  DB  castor  oil  with  28  parts  of  tolylene  diiso- 
cyanate; the  mixture  is  cured  12  to  18  hours  at  158°F.  The  specific  gravity 
is  1.03.  The  complex  shear  modulus  as  a function  of  temperature  and  pressure 
is  given  in  reference  8.  This  polyurethane  can  also  be  classified  as  a stiff- 
lossy  polymer  when  compared  with  the  low  shear  modulus  rubber  compound.  As 
seen  from  figure  8,  the  additional  stiffness  of  the  elastomer  significantly 
raises  the  fundamental  array  resonant  frequency,  and  the  damping  due  to  the 
elastomer  decreases  the  Q value  of  the  array  resonance  or  broadens  the  insertion 
loss  resonance. 


CONCLUSIONS 


Good  agreement  has  been  found  between  the  calculations  of  the  analytical 
model  and  measurements  of  insertion  loss  from  tightly  packed  compliant  tubes 
in  viscoelastic  coatings  in  which  the  major  axes  of  the  compliant  tubes  are 
perpendicular  to  the  direction  of  the  normally  incident  wave.  High  shear 
modulus  polymers  tend  to  both  raise  the  frequencies  and  decrease  the  quality 
factor  of  the  array  resonances.  Near  the  array  resonance  frequencies  of  the 
tubes  in  a low  shear  modulus  material,  the  resulting  low  impedance  of  these 
coatings  tends  to  reflect  much  of  the  incident  energy.  With  increased  stiff- 
ness of  the  encapsulant,  more  energy  is  transmitted  through  the  layer  at  the 
array  resonance.  In  general,  the  effects  of  a high  shear  modulus  polymer  are 
temperature  sensitive. 
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Figure  6.  Comparison  of  Calculation  and  Measurement  of  Insertion 
Loss  for  0.375  in.  (9.52  mm)  Lexan  Plastic  Tubes  in  a 
0.45  in.  (11.50  mm)  Spaced  Grating  Encapsulated  with  a 
0.5  in.  (12.7  mm)  Layer  of  Polyurethane 
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Further  improvements  in  the  model  should  include  generalizations  of  tube 
orientation  and  angle  of  the  incident  radiation,  plus  the  inclusion  of  the 
analysis  of  staggered-tuned  and  multiple- layered  compliant  tube  arrays. 

These  more  complex  arrays  increase  both  the  magnitude  and  bandwidth  of 
insertion  loss.18  The  implications  of  the  complex  elastic  moduli  treat- 
ment have  been  recently  investigated  by  Gaunaurd.19  He  showed  that  this 
model  is  equivalent  to  the  Kelvin  model  with  harmonic  time  dependence 
assumed.  In  future  work,  comparison  with  other  viscoelastic  models,  such 
as  the  Maxwell  model,  should  be  given  consideration. 
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